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ONS  AKE  THE  DOMINANT  EIGENVALUE  PR05I.EM 


JEFFERY  J.  LEADER 
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dominant  eigenvalue  problem  are  analyred  from  a  ci serene 
dynamical  sysvems  perspective.  It.  is  shown  that  the  method 
can  ext.ract  more  information  than  the  standard  power  method 
but  at  greater  computational  cost. 
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Trie  pcw&r  iteTCLtioTi  is  the  matrix  i'Leratiori 

V  =  B^V  .-’liV  ii  Ci  .  !> 

n-*-l  n  ri 
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to  the  pcj'u>€^r  method  f  op  f  i  Pidi  Pig  the  dominant  eigePi  value  of  a 

peal  mio  p  i  x  % 

\  =  Ax  ,  i;  -  V  yp  C 1  .  E: J 

"  ri-r  J  7*)  rt-4- 1  n- 1  r, 

vhepe  A  is  a  real  maipix  with  a  dominant  eiaen value*  v  is  an 
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kr-  li:  i  ITlc 


okr"pi t  c jp  associ  «=!. t e ' 


domi  rian*  eiaenvalue  of  A,  arid  p  is  an  element 

r 

t  pp  Ol  -ep  t  y  t^hat. 

rni  "  r'r..iiia. 


With  the 
'v  wi  t  h 


C  s  l1  , 


ei  aePiVa  j  uev.  s  j  oi 


’We  will  show  that  although  the  power 
L  .  j  V  j  o e Pi e p  cT' j.  j.  Sx  owe  r  thapi  t-htr  pO'^t^j  method 

a  I  i  pro  Vi  de  ex>»pa  irifopm<nt  iopj  ci*.joUs.-  tht;  q  c^irii  pi  c*  n  t 
a  rria  trix  in  certain  cases.  Ve  take  a 


1  ^  cx^  d  t  e  V-  ^  Li  V  j  i  o.  u  t  X  cj,  X 


Q tr-  e  vpii— t  P  ,  o  P^pr  L.-a  L  Ti  ,  Vx  e  v«.'.L  ' 

swsoerPi  dPiO  i  Pipui  p  i  PiO  as  to*  ohe  Piauur  e  ol  its  i  i  iTix  o  sets 
C  at.tr  act  or  in  various  cases  Cin  the  spirit  of  [9]J. 

Thie  iteration  Cl.  ID  is  considered  in  a  different,  context  as 
a  special  case  of  the  CR^->[R^  map 


V  =  AxV  -I-  ExV  XlIV  II  Cl  .  3D 

n+l  n  n  n 

in  [101;  based  on  work  in  [El  Calso  reported  in  [31D.  Further 
deta  i  1  s  opi  the  i  ter  at  i  opi  C  1 . 3D  may  be  f  oupid  in  [  1  C-> ,  i  E  ,  1  5  j,  and 
the  forthcomiing  [41.  Although  C 1 . 3D  is  opil  y  a  linear 

perturbation  of  the  vvel  1 --behaved  iteration  Cl.  ID,  it  exhibits 


d 


course  ,  C  1  .  1 

can 

al  so 

be 

VI  ewed 

as  Cl. 

,  3j  wi  th  a  change  of 

nor  mal i zati on , 

and 

much 

i  s 

k  nown 

about 

the  numerical  method 

gi  ven  b  ^’  C  j.  . 

C  see 

al  so 

C8, 

p.  353, 

j  j  . 

a. 

The 

Power 

Method 

1  l"ii0  Cl.t—.,'  l~-»ciS  xl  A  is  c*. 

nondefeciiv^  matrix  with  a  dominant  eioenvalue,  say  >.  ,  and  ^> 

'  1  o 

has  a  nonzaro  projection  on  an  ei  convector  associated  wdth 
this  don.inani  eigenvalue,  then 

P  — >  /v  as  n  — ^  cXi 

n  1 

and  V  converges  to  an  eiaen vector  associated  with  X  and  'with 

n  1 

unit  /  nor  mi.  If  v  does  not  have  a  nonzero  component  along 
a>  o 

an  eigc  -  n  vec  "Lor  associated  'with  X  and  infinite  precision 
arithme"  ic  is  used  then  we  must  consider  the  eigenvalue  of 
largest  modulus  along  w'hi  ch  does  have  a  nonzero  component. 

In  actual  computations,  how^ever  ,  a  comiponent  along  an 
eigenvector  associated  wdth  would  almost  certainly  be 

introdij-ed  eventually  and  magnified  in  successive  iterations 
[1,  p.l45].  The  same  results  are  found  if  A  is  defective 
Cconsidering  now'  principal  vectors  [7,  p.  31  rather  than  just 
ei  genvector s J  but  the  convergence  is  much  slow'er. 

When  the  dominant  eigenvalue  is  real,  the  method  converges 
to  a  fixed  point.  However,  w'hen  the  dominant  eigenvalue  is  a 
complex  conjugate  pair,  the  method  generally  fails  to 
converge.  Methods  exist  to  recover  informiation  in  such  cases 


[6,  p.co7]  but  they  tend  to  be  somew'hat  involved. 


tTiiCT-ji  V  «rjjfer  icr-<^-L,  ‘w-lj 

O 

./  —  A-^  V  "  '  '•' y  ■. 

n-t*  1  n  n  00 

and  :  n  1  ni  s  f  ormul  at  i  on  nf^a-d 

r» 

actual Iv  needed  Cto  estimate  A  ‘ 

1 

IS  the  same  as  Cl. 11  except  for 
the  nor Jif-il  1  zati  on .  For  this  re 
quantity  ^  in  C 1  .  SI  as  the  si^r 

n 


.^.4“  \ - - 

o.  v^ii<nii^t=;  »vr-  I  *o  V  er 


V  —  A^y  >■■'  ti  A^  v  li 

n-^l  n-^-l  ‘  n  00 


not  be  calculated  until  it  is 


Then  the  iteration  for  v 

n  +  l 


the  particular  /  norm  used  in 

r 

ason  we  sometimes  refer  to  the 


ted  /  norm, 
a 


Conic  C'rbits 


If  h  IS  nonsinciular  then  the  points  V  ,  V  ,  .  of  the 

1  ?  3 

orbit  of  Cl.il  all  lie  on  a  conic  defined  by  the  matrix 

G  = 


1 c .  a 1 i  1 .  For. 


V  =  1 


V 

n-t-i 


C  V  OxGxCBxV  /IIV  ID 

n  r»  nr. 

v'^xCB^xGxB:)xV  /IIV  II 


v’^xjxv  /liV  |!‘ 


for  every  n>l  and  for  any  V  which  is  nonzero.  Clearly  G  is 

^  o 

positive  definite  symmetric,  and  so  the  points  •  • 

must  all  lie  on  the  hyper ell  ipse  defined  by 


V  =  1 


C3.  11 


1  n  [R  . 


If  B  IS  singular,  a  similar  result  holds.  In  order  to 

handle  Limul  taneousj.  y'’  both  the  ca-se  v'here  B  is  simple  and  the 
case  where  B  is  defective  we  state  the  result  in  terms  of 


.<1 


u’  .  W.  i  V. 4^  .-.1  1  .......  »-■.  --  XV'..- - .  ♦-  .-.  .- 

iicr  i  i<ji  ^  ^  cr^ji. 


Theorem  1:  Suppose  B  is  a  real  square  mxm  matrix  wi-Lh  0<q<rr! 

nLil  1  eigenvalues.  If  q>  O  and  has  a  nonzero  component  along 

a  principal  vector  associated  with  a  nonzero  eigenvalue,  then 
all  orb  ts  of  the  power  iteration  Cl.lJ  are  constrained  to  a 
hyper  el  1  i  pse  in  C  m-q*^  -di  mensi  ons  C  f  or  al  I  but  f  i  ni  tel  y  many 
n!)  .  Oth  erwise,  the  orbit  reaches  the  origin  in  finitely  many 
.ti  cns . 

Proof:  First,  not.e  that 

V  =  B^V  /ItV  II 

1  o  o 

V  =  P.X  V  /•  j(  V  II 

2  11 

^  B^CB^V  --liV  xjiV  li  It 

o  o  o  o 

=  /liB^V  II 


aild.  in  gene  raj. 


for  n>l 


Nov.'  let 


V  =  /IIB'’'"'^xV  I 

n  U  O 


C  3.  BD 


j  =  R 

be  the  Jordan  normal  form  of  B  for  some  nonsingular  R. 
Substituting  this  into  Cc>.S?j  gives 

V  =  /iib’’''*»<v  II 

n  O  O 

o  o 


C3.  31) 


where  Z  =R  and 

o  o 


a  =  1/IIB'^  II 

n  O 

is  a  scalar  Cfor  each  nD .  Cl  earl  v,  if  Z  has  no  nonzero 

o 

component  along  a  principal  vector  of  J  C equi val entl y ,  if 

has  no  nonzero  component  along  a  principal  vector  of  EJ  that 

IS  associated  with  a  non-nuLl  eigenvalue,  the  term  J  '^Z  in 

o 


o 


”v  •c'  c  t/-  o  r  . 


r>-*- 1 

und^f  1  r-is  ci  a^nd  -Lhv^  i -t^^r  <=*  1 1  o)  i  s-Lops.  Othns^r  wi  ,  for  n 

sufficiently  large  Cit  suffices  that  ,  all  Jordan  blocks 

in  J  associaced  with  a  null  eigenvalue  will  have  become  blocks 
of  enti’^ely  zeros  in  J^’  Csince  these  Jordan  blocks  are 
nilpctentD.  Now,  the  pri nci pal  vectors  belonging  to  a  given 
Jordan  block  do  nol  interact  with  the  remaining  principal 
vectors,  in  the  sense  that  if  u  is  a  principal  vector 


j 

associated  with  J  C  ,  a  Jordan  block  of  J,  then  J^’^v 


^  J 

involve!  only  a  linear  combination  of  principal  vectors  of  J 
uhat.  ar»  also  associatf^d  with  J  .  Therefore  for  n  large 

V 

enough  ’hat  a.ll  nilpotent  Jordan  blocks  have  become  entirely 
zero  submatrices,  the  vector 


o 

can  written  in  terms  of  a  basis  consisting  of  only  the 

rernainirg  i  m-ql)  principal  vectors.  Thus  the  iteration  lies  in 
a  C  m-q") -di  mensi  onal  subspace  of  [R"^,  and  a  suitable  change  of 
variables  car  then  be  used  to  transform  the  iteration  into  one 
of  the  f  or  iTi 


V  =  C^V’  kliW  II 

n-^  1  n  n 

where  V.‘  is  a  Cm-qJ -vector  for  every  n  and  C  is  a  real 
nonsi  ngul  ar  Cm-q!)xCm-qJ  matrix.  Hence,  in  this  subspace,  the 
iteration  is  constrained  to  the  hyperellipse  determined  by  the 
rriatr  i  x 

T  -1 

CC^C  J 

Cas  w'as  shown  previously  for  the  nonsingular  caseJ  for  all  but 
finitel;^^  many  n.  d 


We  em^phasize  that  this  is  not  an  asymptotic  result;  after  a 


hyper  G-l  I  1  pse-  C  assuming  infini*Le  precision^.  We  now  wish  1.0 
look  at  the  orbits  on  the  attract.!  ng  hyper  el  1  i  pses . 


4.  Limit  Urbits  for  Nondefectiv'e  Matrices 


Suppose  that  B  is  nondefective  and  nonsingular,  and  le: 

be  a  set  of  eiaenvectors  associated  with  the 


1 


eioenvalues  A.  .  .  .  .  ,  >  ,  respectively,  with 

1  rn 


xj  >  I  xj  > 


A. 


any  V  elk.  may  be  written 
o 


V  =  y  a  X 
o  ^  \  \ 


vAiere  a  is  the  component  of  along  x  .  Computing  the  power 


iteration  C  1  .  1 1)  by  C3.2D  gives 

m 

V  ‘ 


[m  -  TTi 

r  Oi  br  I  r  CV  x'^'”*v 

^  I  V  V  I  ^  V  V  I 

V  =  J  V  =  1 


C4.  1  J 


or  TiiJLj.  . 


Now  suprjos-  vhat  B  has  a  C  repeated!)  domi  nant  eigenvalue  a 
with  multiplicity  l<r<mi,  that  is, 

X  =  X  =  •  •  •  =  X 

12  r 

and 


IX  I  >  |X  I  > 

1  r-^1 


>  IX 


Removi  no  a  factor  of  from  the  numerator  and  X^  ^  from  the 

1  1 


denominator  of  C  4  .  1 !)  gives 

V  =  x"  -  fo  X  + 
n  1  1  1 

.n-l  P 

X  *  j  a  X  •  • 

1  It  1  1 

and  cl  er-r  ly,  in  the  limit 


c]  * 


+  a  X  €>C  X  xX 
r  r  r+i 


+  a  X  +  OC  X  xX  D  t! 

r  r  T  +  l  1  ll 


V  — X 
n  n 

wher  e  X  i  s  a  uni  t  vector  i  n  span<  x  ,  .  .  .  ,  x  >  ,  and 

1  r 


7 


J.  =  X  -ex  /  |X  n'"'  C4. 

n  111 

Cprcvidei  that  at  least  one  of  <:>  ,.  .  .  ,o  is  nonzero; 

1  1 

otherwise,  we  begin  the  analysis  anew  bv  consider  inc  X  D. 

r-H 

Hence  i-  X  >0  the  iteration  tends  to  a  fixed  point,  and  if 

1 

X^<0  the  iteration  tends  to  a  symmetric  Cin  the  origin^  period 
two  cycle  on  the  points  ±X  x. 

Now  suppose  that  B  has  a  real  positive  eigenvalue  X^  of 
multiplicity  p,  so  that 


X  =  X  = 
1  ? 


=  A 


and  a  real  negative  eigenvalue  -X^  of  multiplicity  q,  so  that 

=  X 

pi2  p-*-q 

and  suppose  further  that 

■  •  ■  s  I A  I 

m 

P^r  ooeeQi  ng  as  bel  or  e  >  'we  have  that 

<DC  X 


X  =  X 

p-^l  P-+2 


i  A.  1  /  I  X 

1  '  p-^q-*-i 


^  F  P  q  *N 

=  >'  •  r  Oi  >  ±  T  a  >  +  ©c  X 

r,  I  V  I  ^  XI  p-^q-I  1  J 


P 

Yi  o  y 


C  wj  t,  h  / 


V  1 

1  =  1 

oi  ven  in  C  4 


V  =  p  1 
p+  q 

E 

V  i 


OC  X  -AX  D 

p-^q-^l  1 


'  J  Vv'hfc^'r  the  d  IS  po>Sl  tl  Vfc=-  W’hen  n  IS 

r  I 

even  and  negative  when  n  is  odd  in  the  numer a uor ,  and 
contrariiy  in  the  denominator.  In  the  limit,  we  have  that 

C  approx:  matelyj 

V  =  y  -Qy  +  yjyWy  -  y  II 

n  n  1  2  1  2 

V  =  y  -Cv  -  y  D/liy  +  y  II 

n-l  ‘  'i  "^2  1  2 


wner  e 


y  =  Y  &  X 

1  ^  V  t 

X  =  1 

p-*-  q 

y  “  X]  ^ 

2  t  I 

V  =  p  1 

Since  X  13  positive,  y  =X  for  all  n,  and  so  we  have  a  period 

1  n  1 


o 


V  =X-Cy  -i-  y  J  y\\y  -vl 

V  =  X  -Cv  -  y  D/lIv  +  y  II 
1  *^1  2  ''1  2 


which  in  general  is  not  a  syr\jretr  ic  orbit  C  i  n  the  on  gi  n  J  .  We 
will  not  in  general  have  11 V  II  =X  in  this  case,  but  note  that 

n  1 


!iV  II  =  |X  I  •  liy 

ri  11 


y  li  X  li  V  -  y  I 

'^Z  '"1  2 


iiv^^ji  =  ■"  ^2'' 


so  that 


li  V  II  •  li  V  ,i  =  I  X  I  ^ 

n  1  1 

1 n  the  1 1 n  t  of  1 ar  ge  n ;  the  modul us  of  the 
ei  genvo.l  -le  i  s  the  geometr  i  c  rnean  ol  the  norms  ot  owo 
success!  ve  i  ter  ates  C  j  n  the  1  i  m:  t'j  ,  Note  that  the  asymmetry 
of  the  orbit  allows  us  to  distinguish  this  case  f rom  the  case 

ax  ways  gi ves 

_  r  01  i- ^  jTtv^sf  U‘i  th^  ti  m*zr'  C 1  .  fer .  wht^n  tht? 

esul t: r g  asymptotic  orbit  is  indeed  asymmetric!. 

Now'  1  t  us  suppose  that  B  has  a  complex  conjugate  pair  of 


W'htrr  tr  A  i 

1 


a  ricrgaui  v^  doriii  ^iant  ti  gtrri'v''ai  Uc-'  c  ’whi 


sym 


I  jl  I  mil  'win' 


eiaenvalues  that  is  dominant,  i  .  e.  that  A  and  \  are  comiolex 
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u  on  J  u  wg a  t  o  s 


anc 


A 


Fr  w^mt  C  4- .  1 !  w'^r-  havfc- 
V 


,  ,  f  -V 

=  X  '  ae  X  +  oe 

±11  1  2 


ri&  i.  V  ^ 

X  +  OCX  xX  ! 

2  a  1 


v<r»-i>@  -v(n-l>@  N  v  ^  ri-1 

ae  x+oe  X  +  OCX  xX  ! 

1  12  2  3  1 


whiere  b'-’argCX  !.  In  the  limit,  this  becomes 


V  = 

n 


,  .  ,  f  x.’nO  -vn0 

=  X  *  ae  X  +  o(e  x  e 

‘  1  ‘  1  1  2  z) 

Un-i)^  -i,(n-l>6' 

o  e  X  +  a  e  x 

1  12  2 


C4.  3! 


Thus  if  ©is  such  that  expC  v©J  =e>qDC  vC  n+1 !  6!  for  some  n  Ci.e. 


if  X^  and  X^  are  both  in  for  some  n!  then  we  have  an 


the  or  ba  t  is  aperiodic  on  the  underlying  hyp-er  el  I  i  pse .  If  the 
orbit  is  period  n  then  w'e  note  that  the  geometric  mean  of  the 
norms  o:  n  consecutive  iterates  tends  to  \\  |,  as  is  easily 
seen  by  writing  ouf  tne  proouct  oi  n  iterates  and  noting  the 
ca  ncel  1  a  t i on . 

In  a  similar  w'ay,  if  E  has  tw'o  pairs  of  complex  conjugate 
eigenvalues  of  equal  modulus  and  all  other  eigenvalues  of  E 
have  lesser  moduli,  then  w'e  have  a  situation  like  that  of 
C4.3’i  save  that  there  is  an  additional  angle  to  be  considered. 
HefiCf-  if  the  f  i  rsT  pair  alone  w'oui  d  give  a  pe*' iodic  orbit  of 
period  s,  and  the  second  pair  alone  ’would  give  a  periodic 
orbit  of  period  q,  then  the  orbit  of  the  iteration  Vr'ill  be 
periodic  with  period  n=lcmCp,ql).  If  either  pair  alone  w'oul  d 
give  an  aperiodic  crbit,  then  the  orbit  is  aperiodic.  Again 
the  igeoretric  mean  of  li  consecutive  iterates  tends  to  I  I 
wh'^n  t,h'-  orbit  is  asymptotically  period  n.  The  obvious 

geri*=r  al  1 ‘-an  on  holds  for  more  than  t’wo  complex  conjugate  p^airs 
of  equal  moduli.  In  particular,  in  an  even-dimensional  space, 
say  of  dimension  Er  ,  choosing  B  to  be  a  Sr  x2r  matrix  with  r 
complex  conjugate  pairs  of  eigenvalues,  all  of  equal  miodulus 
and  such  that 

X  ,'"'  e  IR 

V 

for  all  i =1  , .  .  .  , Sr  and  n=l  , S, 3, .  .  .  gi ves  a  method  for 
generating  a  sequence  of  points  on  the  ellipse  given  by  C3.1D 
by  iterating  dl.lJ  with  an  arbitrary  nonzero  V^. 

If  in  addition  to  some  number  of  complex  conjugate  pairs  of 
eigenvalues  wuth  equal  modulus  there  are  some  number  of 


i  n 


t-',  “  A  ^ 

V_y  W  J.  X  '*  5^  ^ 


.•— f  ,  (•w,  |.  «w.  ^  t*  -  1  ^  *-^i  f^i  ♦  »  T  ♦  ♦  ~— *. »—  +  A  W,  *'^.  ^  »•  /-'.  .•— ,  .* J  *  «  1  «  •  ^- 


then  it  is  immediate  that  the  iteration  is  aperiodic  if  any 
compl  e>:  conjuaate  eigenvalue  gives  an  aperiodic  orbit,  and 
periodic  otherwise,  w'ith  p^er  i  od  given  by  the  least  corTLinon 
multiple  of  the  individual  complex  conjugate  pairs  and  the 
period  tv.^o  due  to  the  negative  real  eigenvalues,  if  present. 

We  have  established  the  f ol  1  owi no  theorem: 


1 HEOREW 


r-  :  m  ■;  H  &  r 


V  —  El"^  'v'^  /■  li  \  (! 

n-*- 1  n  n 

vbere  &  is  a  real  mxm  nondefective  matrix  and  V  is  a  Given 

o 

nonzero  m-vector.  Suppose  that  has  a  nonzero  comp*oncnt 

along  an  eigenvector  of  B  which  is  associated  with  an 

eigenvalue  of  B  w'lth  maximum  modulus.  Thien  all  iterates  for 
n>rn  are  constrained  to  an  ellipse  in  a  subspace  of  [? "’  of 
dimension  rankCBD,  and  the  asymptotic  behaviour  of  the 
1 1 er  ates  is  as  f  ol 1 ows  : 

i'J  if  B  has  a  rriultipie  real  dominant  eigenvalue,  then  the 

iteration  tends  to  a  fixed  p>oi  nt ; 

tO  if  B  has  a  multiple  negative  real  eigenvalue,  then  the 

iteration  tends  to  a  symmetric  period  tw^o  cycle; 

iiiJ  if  B  has  both  positive  and  negative  real  dominant 

eigenvalues,  then  the  iteration  tends  to  a  period  two  cycle; 

tvD  if  E  has  a  complex  conjugate  dominant  eigenvalue,  then 
the  iteration  tends  to  a  period  n  cycle  if  the  eigenvalues  are 
a  multiple  of  an  ntdi  root  of  unity  and  is  aperiodic  otherwise; 

v3  if  B  has  multiple  complex  conjugate  dominant  eigenvalues 
and  multiple  real  dominant  eigenvalues,  then  the  iteration 
tends  to  a  periodic  orbit  with  period  equal  to  the  least 
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eigenvalues  when  all  eigenvalues  give  rise  to  periodic  orbits, 
and  IS  aperiodic  otherwise. 

5.  Limit  Orbits  for  Defective  Matrices 


TTie  ^  oove  results  are  essentially  unchanged  if  B  is 
defective.  Let  B  be  defective  and  let  x  be  an  eigenvector 
associated  wi  t  hi  a  dominant  eigenvalue  X.  Let  y  be  a  linear 
combination  of  principal  vectors  of  B  associated  with  X.  Then 


the  iter  atioCi 


^  y 

converges  to  x  as  OC 1  xnl-  TIG,  p*.  SS]  ,  Ci6,  p.  SSo]  .  Since  the 
pr 3 nci p^ 1  vectors  corresponding  to  distinct  eigenvalues  are 
noni  nter  ac  ti  ng ,  it  is  clear  that  the  qualitative  results  of 
Theorem  E  are  unchanged.  although  the  convercenc.e  to  the 
asymptotic  orbits  may  be  exceedingly  slow'.  Hence  TTieorem  E 
remains  valid  if  w'e  allov/  E  to  be  defective,  and  require  that 
have  a  nonzero  projection  on  a  principal  vector  of  E 
associated  w'l  th  an  eigenvalue  of  maximum  modulus.  If 

J=R  IS  the  Jordan  normal  form  of  B  with  the  dominant 

eigenvalue  C  of  multiplicity  p!)  in  the  first  blockCsJ  along  the 
diagonal ,  then  these  vectors  have  the  form 

L 

C  i  =1  ,  .  .  .  ,  pJ  ,  where  <e  is  the  natural  basi  s  f  or 

\  I  =1 


The  Power  Iteration 


coiij  LiTisi  C'l  or 


j.  i  tr ji 


met  hod  for  the  numerical  del  er  mi  nahi  on  of  t.he  modulus  of  the 
doiTii  nani.  eigenvalue  of  a  real  matrix  when  the  resulting  orbit 
is  periodic.  If  the  orbit  has  period  n^l ,  then  the  modulus  is 
approxi  iTiatel  y  the  geometric  mean  of  the  Euclidean  norms  of  n 
consecutive  iterates.  Addi t i onal 1 y >  the  argument  is  such  that 
the  power  of  both  the  eigenvalue  and  it’s  conjugate  are 

positive  real  numbers,  so  that  the  desired  eigenvalue  K  is 
gi ven  by 

X  =  |\  I  CO  C6.  ID 

V 

f  or  somt  1  ~1  ,  .  .  .  ,  n ,  where  co  ,  .  .  .  ,  co  are  the  nth  roots  of 

1  n 

unity.  If  the  orbit  is  aperiodic,  the  dominant  eigenvalues 
are  complex  conjugates  and  fail  to  satisfy  C6.1D. 

As  an  example  of  the  use  of  the  powder  iteration,  consider 


tfi»rT  mesitri 


■  (ni 


C6.  £D 


which  s  eigenvalues  l/d.  In  C14  p.  9S]  it  is  shown  that  the 
pc-iwer  method  C 1  .  ED  app»l  i  ed  to  this  iteration  C  wi  th  initial 

T 

vector  (1,1D  D  settles  into  a  period  two  cycle  on  the  two 
vector  s 


^  =  [-a 


i  .  e.  th«:  method  fails.  Thie  pow’er  iteration  Cl.  ID  applied  to 
this  matrix  yields  asymptotically  the  period  two  cycle  Cfrom 
the  same  initial  vector D 


-,■-(3 


X  =  C  S/Vo? 
2 


[i] 


Since  this  is  an  asymmetric  period  two  cycle,  from  Theorem  2 
it  foil  CVS  that  there  are  both  positive  and  negative  dominant 
eigenvalues,  and  that  the  modulus  of  these  eigenvalues  is 


I 


I  A  I  =  C  llx  II  •  ll>  lO^"  ' 

1  2 

^  [  C  ^£0  •  C  2-/£/-/5:j  ] 

-  Vci 

Hence  the  eigenvalues  of  A  are  ±V£ ,  as  expected.  Si  itia  1  ar 
results  would  be  obtained  using  any  norm  in  place  of  the 
Euclidean  norm  in  Cl.lj. 

7.  Conclusions 

The  p  'wc^r  iteration  is  certainly  slower  than  the  power 

Jitethoc  Cdue  tcj  the  need  to  calculate  a  Euclidean  vector  nor  ni 
rather  '  han  simply  locating  an  element  of  the  vector  with 
jTiaxi  mum  modulus^,  and  additional  information  can  be  gained  in 
only  a  restricted  set  of  cases  C  primarily  w'hen  the  dominant 
eigenvalues  are  complex  conjugates  and  real  multiples  of  a 
r  oc'i.  c*f  unityt.  Nonetheless,  in  tnese  cases  it  does  pr  ov'i  de 
useful  information  about  the  eigenvalues,  and  in  more  genera} 
cases  1  may  provide  som?^  insight  as  well.  For  example, 

inspecti''ni  of  the  elliptical  orbits  C  Theorem  ID  can  be  used  to 
provide  information  about  the  existence  and  multiplicity  of 
null  ei  * 'enval  ues .  For  these  reasons  the  power  iteration  may 

be  useful  in  certain  circumstances.  Of  course,  the  power 
i  ter  at  i  c  nXpower  method  is  in  some  sense  the  basis  of  most 

iterati\e  met. hods  for  the  eigenvalue  problem  [15]  and  so  this 
analysis  may  be  useful  in  the  analysis  of  more  practical 
algorithms  for  this  problem. 

We  mer.iiori  t-hat  computergr  apihi  caJ  evidence  seemis  to  indicate 


14 


that  the  unsigned  power  method 


^  !i  ,  .  1,  r  T'  -1  " 

V  “  '  "  >  "  ^  »  .  a. 

r>-H  r>  r>  00 

hos  the  orooerty  that  all  iterates  lie  on  a  /  conj  r-  C  wi  t  h, 
resp^ect  to  soiiie  rotation  of  the  axesD  asymptotically.  As 
noted  in  [8.  p.  368]  ,  however,  the  use  of  the  Euclidean  norm  in 
C 1  .  1j  greatly  facilitates  the  analysis  of  the  power  iteration 
Csee  also  [5,  p.  351  3D  and  we  have  been  unable  to  show  a 
cor r esponda ng  result  for  the  unsigned  power  method  C7.1D.  The 
analysis  of  cases  in  Theorem  S  did  not  depend  on  w'hich  norm 
w'as  usel  in  Cl.  ID  and  as  such  it  holds  for  C7.1D  as  well; 
therefore  the  comments  in  £?6 .  exceolina  those  about  the  Z 


conic  orbits,  are  equally  appiicaL^i 


.ne  unsignec  pjower 


method,  which  is  no  more  costly  than  the  usual  power  method 
C  requiring  only  that,  an  extra  absolute  value  be  taken — but 
this  would  be  done  during  the  cor  r  espondi  nq  searck,  for  p  D. 

r> 

A]  though  usi  no  C7.  ID  i  ri  pi  ace  of  C 1 . 2D  means  that  the 
algoijthr  will  not  co^. verge  Cin  th  -  usual  senseD  in  the  case 
cT  a  dr-.-  inant  negaut  ve  real  eigenvalue,  the  method,  p>r  oper  1  y 
interpr^T  .  ed  Csee  also  [6.  p.257jD,  would  avoid  some  of  the 
pr  obi  em '  enr  our  u.  er  e-h  wher-i  using  the  power  method.  For  the 

matrix  A  given  in  C6.2D,  the  unsigned  power  method  C7.1D,  w'l  th 

T 

initial  vt-ctor  Cl  ,1D  ,  gives  a  period  two  cycle  on  the  vectors 


■■ .  a 


■  C) 


and  aor^i  n  the  fact  that  this  is  an  asymmetric  period  two  orbit 
implies  that  the  eigenvalues  of  A  are  ±  j  X.  |  ,  where 

|X  I  =  Cllx  II  •  llx  II 

1  CO  2  00 

=  C2-1D 

=  VEa 

as  exp-rcted.  This  type  of  reasoning  could  easily  be 
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